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ABSTRACT 

Beyond reionization epoch cosmic hydrogen is neutral and can be directly observed 
through its 21 cm line signal. If dark matter (DM) decays or annihilates the corre- 
sponding energy input affects the hydrogen kinetic temperature and ionized fraction, 
and contributes to the Lya background. The changes induced by these processes on 
the 21 cm signal can then be used to constrain the proposed DM candidates, among 
which we select the three most popular ones: (i) 25-keV decaying sterile neutrinos, 
(ii) 10-MeV decaying light dark matter (LDM) and (iii) 10-MeV annihilating LDM. 
Although we find that the DM effects are considerably smaller than found by previous 
studies (due to a more physical description of the energy transfer from DM to the gas), 
we conclude that combined observations of the 21 cm background and of its gradient 
should be able to put constrains at least on LDM candidates. In fact, LDM decays (an- 
nihilations) induce differential brightness temperature variations with respect to the 
non decaying/annihilating DM case up to AST}, = 8 (22) mK at about 50 (15) MHz. 
In principle this signal could be detected both by current single dish radio telescopes 
and future facilities as LOFAR; however, this assumes that ionospheric, interference 
and foreground issues can be properly taken care of. 

Key words: intergalactic medium - cosmology: theory - diffuse radiation - dark 
matter 



1 INTRODUCTION 

Cosmic recombination at z ~ 1000 left the gas in the Uni- 
verse in a nearly uniform, dark, neutral state which we cur- 
, rently denote as the Dark Ages of the Universe. Gas re- 
' mained neutral until the first luminous objects emerged at 
z ~ 20—30, leading to the complete reionization of the in- 
tergalactic medium (IGM) at a later, yet unknown, epoch. 

The investigation of the Dark Ages is one of the fron- 
tiers of modern cosmology and a new generation of radio 
interferometers is currently in development to study the red- 
shifted 21 cm hyperfine triplet-singlet level transition of the 
ground state of neutral hydrogen. Instruments such as the 
LOw Frequency ARray (LOFAR), the 21 Centimeter Array 
(21CMA), the Mileura Wide-field Array (MWA) and the 
Square Kilometer Array (SKA), are expected to reach the 
sensitivity required to map the HI distribution at angular 
resolution of the order of a few arcminutes (e.g. Pen, Wu, & 
Peterson 2004; Bowman, Morales, & Hewitt 2005; Kassim 
et al. 2004; Wyithe, Loeb, & Barnes 2005). 

The main scientific goal of such instruments is to per- 
form an accurate tomography of matter during these remote 
cosmic epochs, and reconstruct in detail the latest phases 
of the reionization history. However, this might not be the 
only major achievement of such experiments. For example. 



an analysis of the features of HI 21 cm maps from the Dark 
Ages could also allow to put constrains on the existence and 
nature of decaying and annihilating Dark Matter (DM) . In 
fact, if DM actually decays or annihilates, the injection of 
high energy photons into the IGM would heat and ionize the 
neutral gas, leaving an imprint on the 21 cm background sig- 
nal which could be directly observed (Shchekinov & Vasiliev 
2006; Furlanetto, Oh, & Pierpaofi 2006). 

The importance of DM decays and annihilations on the 
reionization history has been analyzed in detail by several 
authors. Most authors have come to the conclusion that the 
effects of DM on the overall reionization process are rela- 
tively small and they cannot compete with those induced 
by more conventional ionizing sources as stars and (mini-) 
quasars. However, as we will explain in detail later on, ra- 
diation from decaying/annihilating DM might change the 
thermal and ionization history of the gas to such an extent 
that differences with the case in which DM is not considered 
are large enough to produce a clear signature in observable 
21 cm signal. Assessing the amplitude of these deviations 
from the standard scenario requires a careful study of com- 
plex physical processes, and in particular of the amount of 
energy that can be transferred from the decay /annihilation 
products to the IGM. A recent study (Ripamonti, Mapelli, 
& Ferrara 2006a, hereafter RMF06a) finds that only a rela- 
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tively small fraction of the injected energy is effectively ab- 
sorbed by the IGM and goes into heating and reionization. 
As a result, the effects on the 21 cm background signal might 
be smaller than previously predicted (Shchekinov & Vasiliev 
2006; Furlanctto ct al. 2006); hence, the question remains if 
the DM signal can still be observed by next generation radio 
interferometers. If so, this experiment would constitute a su- 
perb tool to distinguish among different DM candidates via 
their decaying/annihilating properties. In this paper we cal- 
culate the effects of DM decays/annihilations on the 21 cm 
background for some of the most popular DM candidates, 
such as decaying or annihilating Light Dark Matter (LDM) 
and decaying sterile neutrinos. 

The rest of the paper is organized as follows. In Sec. 2 we 
briefly introduce the DM candidates and their effect on the 
IGM; in Sec. 3 we provide the basic equations to study the 
21 cm radiation in presence of decaying/annihilating DM. In 
Sec. 4 we present the results of our calculations, which are 
then discussed in Sec. 5. Throughout the paper we assume, 
in agreement with the 3-yr WMAP data analysis (Spergel ct 
al. 2006), a ACDM cosmology with fi^ = 0.24, = 0.76, 
h = 0.73, ^6 = 0.042, Ho = 100 /ikms"^ Mpc'^ 



2 BASIC PHYSICS 

We are interested in calculating the effects on HI 21 cm 
line signal produced by two among the most popular low- 
mass DM candidates, i.e. LDM (~ 1 — 10 MeV) and sterile 
neutrinos (~ 2 — 50 keV). In the case of sterile neutrinos 
only the decay process is allowed; LDM particles instead 
can both decay and annihilate. We neglect heavier DM can- 
didates (with mass larger than ~ 100 MeV) because previ- 
ous studies (Mapelli, Ferrara & Picrpaoli 2006) have already 
shown that, even assuming that all the energy released fol- 
lowing annihilations is immediately absorbed by the IGM, 
they represent a negligible heating/ionization source for the 
gas. 

Both in the case of sterile neutrinos and of LDM, the 
rate of energy transfer per baryon to the IGM can be written 
as (RMF06a; Ripamonti, Mapelli & Ferrara 2006b, hereafter 
RMF06b): 

Ea:{z) = fabs(z)nDMiz)mDMC^, (1) 

where ttidm is the mass of the DM particle and c is the speed 
of light. The energy absorbed fraction, faba{z), is discussed 
in detail in the following; ti-dm^z) is the decrease rate of the 
number of DM particles per baryon. 

In the case of DM decays, nDM(a) is given by 

nDM(2) ^ , (2) 

TDM 

where noM.o and tdm are the current number of DM par- 
ticles per baryon and the lifetime of DM particles, respec- 
tively. 

For the annihilations: 

hnviiz) ~ i nDM,o-'Vb(0) {cjv){l + zf, (3) 

where A/'b(0) = 2.5 x 10~^ cm~^ is the current baryon num- 
ber density (Spergel et al. 2006), and {av) is the thermally- 
averaged DM annihilation cross section. The values of twdm. 



riDM.o, TDM and {av) depend from the properties of the cho- 
sen DM candidate (see Sees. 2.1 and 2.3). 

2.1 The energy absorbed fraction 

The most delicate term in cq. (1) is represented by fabs(z), 
i.e. the fraction of the DM particle rest mass that is ab- 
sorbed by the gas at a given redshift z; fabs strongly de- 
pends on how the decay/annihilation products interact with 
the IGM. Since accounting for the physical processes which 
govern such interactions is quite complicated, most of pre- 
vious studies assume that: (i) all the energy released by 
DM decays/annihilations is immediately absorbed (Hansen 
& Haiman 2004; Pierpaoli 2004; Biermann & Kusenko 2006; 
Mapelli et al. 2006), (ii) leave fabs as a free parameter (Pad- 
manablian & Finkbcincr 2005; Zhang ct al. 2006), or (iii) 
make a partial treatment of the energy redistribution (Chen 
& Kamionkowski 2004; Mapelli & Ferrara 2005). 

Recently, RMF06a have calculated behaviour of fabs (z) 
in detail, for the most common case in which the de- 
cay/annihilation products are photons, active neutrinos or 
electron-positron pairs. Photons are affected by Compton 
scattering and photo-ionization; for pairs, the relevant pro- 
cesses are inverse Compton scattering, coUisional ioniza- 
tions, and positron annihilations. 

RMF06a found that, if the decay /annihilation products 
are either photons and active neutrinos or pairs, fabs is close 
to the maximum allowed value (equal to 0.5 for sterile neu- 
trinos, due to the active neutrino production, and equal to 
1 for LDM) only at very high redshift {z > 100). At lower 
redshifts, fabs rapidly drops to values < 0.1; also, the higher 
is the mass of the progenitor DM particle, the faster is the 
decrease of fabs- Thus, accounting for the correct fabs{z) de- 
termination dramatically reduces the possible effects of DM 
decays and annihilations on the IGM heating and ionization 
(RMF06a) with respect to previous estimates. 

In this paper, for the first time, we will adopt the correct 
estimate of fabs{z) (as given in RMF06a) in order to eval- 
uate the impact of DM decays and annihilations on 21 cm 
emission. Previous papers (e.g. Furlanetto et al. 2006) have 
assumed a redshift-independent fabs{z), which appears to be 
quite unrealistic (for a discussion, see Sec. 4) and leads to 
optimistic upper limits for the contribution of DM-related 
effects to 21 cm maps. 

Differently from Furlanetto et al. 2006, who do not se- 
lect any specific DM candidate and leave the decay rate as a 
free parameter, we chose to consider three specific DM can- 
didates (sterile neutrinos, decaying LDM and annihilating 
LDM). This choice allows us to give predictions which can 
be more easily related to other DM measurements, such as 
the X-ray constraints on the sterile neutrino mass (Watson 
et al. 2006; Boyarsky et al. 2006a) or the detection of the 
511-keV emission line from the Galactic center (Knodlseder 
et al. 2005). 

2.2 Sterile neutrinos 

Many models of sterile neutrinos have been proposed, with 
mass ranging from eV to TeV. Here, we are interested in ster- 
ile neutrinos as warm DM (WDM) candidates, with masses 
of the order of a few keV (~ 2 — 50 keV). These particles 
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can decay into an active neutrino and a photon (Dolgov 2002 
and references therein). 

The mass (and thus the hfetime) of radiatively decay- 
ing sterile neutrinos can be constrained by the absence of 
any detection of X-ray Uncs consistent with photons duo to 
sterile neutrino decays in galaxy clusters (Abazajian, Fuller 
& Tucker 2001; Abaaajian 2006; Abazajian & Koushiappas 
2006; Boyarsky et al. 2006b) or in galaxies (Watson et al. 
2006 obtained the strongest constraints from the study of the 
Andromeda galaxy). Other constraints come from the com- 
parison between the unresolved X-ray background and the 
expected contribution from sterile neutrino decays (Mapelli 
& Ferrara 2005; Boyarsky ct al. 2006a). 

In this paper we consider the representative case of 
m,/ = 25 keV sterile neutrinos, whose contribution to heat- 
ing is maximum (RMF06a), due to the weakness of the 
available constraints on lifetime for neutrinos of such mass. 
For such particle the upper limits on the lifetime and the 
present number density number per baryon are tdm = 
9.67 X 10^'' s and wdm.o = 1-88 x 10®, respectively (see 
RMF06a, RMF06b). The contribution of other sterile neu- 
trinos masses to the 21 cm line is expected to be comparable 
or smaller than for the case considered here. 



2.3 Light dark matter 

We define as LDM particles all the DM candidates whose 
mass is 1 < rriLDM /MeV < 100. Such particles have been 
suggested as a possible source for the dotoctod 511-kcV ex- 
cess from the Galactic centre (Knodlseder et al. 2005). Ac- 
cording to this scenario, their maximum allowed mass mLDM 
should be 20 MeV, not to overproduce detectable gamma 
rays via internal bremsstrahlung (Beacom, Bell & Bertone 
2004), or even ~ 3 MeV, if we consider also the production of 
gamma rays for inflight annihilations of the positrons (Bea- 
com & Yiiksel 2006). 

LDM can decay or annihilate, producing photons, neu- 
trinos and pairs. We will treat both channels in detail and 
assume that the only decay /annihilation products are pairs, 
an assumption leading to an upper limit in terms of IGM 
heating (RMF06a). 

We consider, as a template, the case of 10-MeV LDM 
particles, which again yields the most efficient heating case 
(see RMF06a). For such particles the upper limits of the 
current number (per baryon), the lifetime and the cross- 
section are noM.o ~ 446, tdm = 4 x 10^® s, and (av) ~ 

2.4 X 10"^^ cm^ s"\ respectively (RMF06a). 



3 EFFECTS ON THE 21 CM RADIATION 
3.1 Lya pumping 

The 21 cm line is associated with the hypcrfinc transition 
between the triplet and the singlet levels of the hydrogen 
ground state. This transition is governed by the spin tem- 
perature, Ts, defined as: 



m 
no 



= 3exp 



(-1) 



(4) 



where no and ni are the number densities of hydrogen atoms 
in the singlet and triplet ground hyperfine levels, and T* = 



0.068 K is the temperature corresponding to the transition 
energy. 

In the presence of the Cosmic Microwave Background 
(CMB) alone, the spin temperature reaches thermal equi- 
librium with TcMB ~ 2.73(1 + z) K on a short time-scale, 
making the HI undetectable in emission or absorption. 

Two mechanisms can decouple Ts from Tcmb : (i) colli- 
sions, which are effective mainly at « > 70 due to the higher 
mean IGM density, and (ii) scattering by Lya photons — the 
so-called Wouthuyscn-Ficld process or Lya pumping (e.g. 
Wouthuysen 1952; Field 1959; Hirata 2005) - which cou- 
ples Ts to the kinetic gas temperature Tk via the mixing of 
the hypcrfinc levels of the HI ground state through interme- 
diate transitions to the excited 2p state. 

The spin temperature is then given by the equation: 

Tcmb + VaTt + ycTk 



Ts = 



(5) 



where Tk is the kinetic temperature, and the Lya and col- 
lisional coupling coefficients are given by the following ex- 
pressions: 

PioT, 



Ha = 



Vc = 



^loTfc ' 
CioT* 
AwTk ' 



where 

Clo = kloUHI + ne7e 



and 

Pio -- 



16 TtJcCTcx 

27 hpUa 



(6) 
(7) 

(8) 
(9) 

^®s~^ is the spon- 



In the above equations ^lo = 2.85 x 10 
taneous emission coefficient of the 21 cm line. Pro is the 
indirect de-excitation rate of the hyperfine structure levels. 

We write Pio = 4Pa/27; in addition. Pa = {4nJacra)/{hpi^a) 
is the outcome of the equation for the Lya scattering rate. 



(10) 



Pa = c di/n(i/) = — / 

J V h J V 

in the case that cr(i^) is a 5 function. A detailed investiga- 
tion of the physics of the Wouthuysen-Field process (Hirata 
2005) finds small corrections to these expressions which we 
neglect here for simplicity. The coefficient Cio in eq. (8) 
is the coUisional de-excitation rate by hydrogen atoms and 
electrons, where k\Q is tabulated in Allison & Dalgarno 
(1969) for different temperatures; 7e, according to Liszt 
(2001), is 

log7e(Tfc) = -9.607 + logTfci/2exp(-(logTfc)*Vl800).(ll) 

for Tfc < 10* K, while 7e(Tfe > lO-^K) = 7e(Tfc = lO^K). 

For the Lya pumping to be effective, a minimum Lya 
background intensity, J a, is required. This is given by the 
condition (Ciardi & Madau 2003): 

-1 TJ -1 -1 

Hz sr , 



J,, > 9 X 10 



-23/1 I \ -2 

(1 + 2) erg cm s 



(12) 



at the redshift of interest. 



3.2 IGM evolution and 21 cm background. 

Next, we want to understand how DM decays/annihilations 
affect the thermal and ionization evolution of the IGM, 
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Figure 1. Left panel: Ts and Tj, as a function of redshift (25-keV sterile neutrino decays, faba ^ ™ RMF06a). Thin solid line: Tcmb', 
thick solid (short-dashed) line: Tg with (without) 25-keV neutrino decays; dotted (long-dashed) line: Tj. with (without) 25-keV neutrino 
decays. Right: Ts and T^. as a function of redshift for 25-keV sterile neutrino decays and for f^bs = 0-5- The lines are the same as in the 
left panel. 



which, in turn, determines the level of the 21 cm back- 
ground signal. The equation that describes the evolution 
of the ionized fraction Xe is the following (see e.g. Chen & 
Kamionkowski 2004): 

dxe _ 1 
~"d7 ~ H{z){l + z) 



(13) 



where Ix ~ Ex/ Ea is the contribution to the ionization rate 
due to DM; Is and Rs are the standard ionization and re- 
combination rates per baryon. Considering eqs. (l)-(3) we 
have that 

(14) 



Ix = Xi{z)^ = fab,{z)Xi{z)^xfx^^^^^ 



for DM decays, and 

E m (? 

Ix = Xi{z)-Er- = fabs{z)x^{z)fx—^nDM,oJ^b{z){av} (15) 

for DM annihilations. In the last two equations Eq = 13.6 eV 
is the hydrogen ionization threshold, nip is the proton mass, 
Xi is the fraction of the energy absorbed by the IGM from 
DM decays/annihilations that goes into ionizations (xi ~ 
(1 - Xe)/3, see ShuU 1979, Shull & van Steenberg 1985), 
fx = ^x{z)/i^b{z), Afb{z) is the baryon number density and 
= 1/tdm. 

The equation regulating the evolution of IGM temper- 
ature can be written as: 



H{z)(l + fH. + x^y 

IXhEx 



(16) 



3fc6H(z)(l + /He-Fa;e) 
where Z-y = {^(TTanTcMB)l^''neC)-: Xh is the fraction of the 



absorbed DM energy deposited into the IGM as heating 
(Shull & van Steenberg 1985) and Jmb = 0.24 is the helium 
fraction by mass. The Lya heating resulting from repeated 
scatterings as the photons are redshifted into the Lyman 
resonances is negligible as it has been recently shown (see 
e.g. Chen & Miralda Escude 2004) . The ionization and tem- 
perature eqs. (13)-(16) are solved using a modified version 
of RECFAST (Seager et al. 1999). 

A third equation is needed in order to compute the 
21 cm backgroundi: the one describing the evolution of the 
Lya background intensity Ja- Following Madau, Meiksin & 
Rees (1997) we write: 



Ja{z) 



Afjihc 
4nH{z) 



eff , , XaExiz) 
XeXpar,2p + XeXHIJeH + -TT-, 



.(17) 



where the first two terms are the contributions from re- 
combinations and coUisional excitations by electron impacts, 
while the third term is the DM contribution. 

In the last equation AfH=0.92J\fb is the number den- 
sity of hydrogen atoms, 0^22 p ^® effective recombination 
coefficient to the 2^P level (Pengelly 1964), which includes 
direct recombinations to the 2^P level and recombinations 
to higher levels, followed by transitions to the 2^P level via 
all possible cascade paths. 

Finally, ^eH ~ 2.2 x 10"* exp (-11.84/r4) cm^s"^ is 
the coUisional excitation rate of HI atoms photons by elec- 
tron impacts; Xa is the net fraction of the absorbed X- 
ray photons from DM decays/annihilations which is con- 
verted into Lya photons (ShuU & van Steenberg 1985), and 
Ti = T/{W^ K). 

Once Ts(z) has been determined through eq. (5), we can 
obtain the 21 cm radiation intensity, which can be expressed 
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Figure 2. Left panel: Ts and Tj. as a function of redshift (10-MeV LDM decays, faba as in RMF06a). Thin solid line: Tcmb', thick solid 
(short-dashed) line: Ts with (without) 10-MeV LDM decays; dotted (long-dashed) line: with (without) 10-MeV LDM decays. Right: 
Ts and Tj. as a function of redshift for 10-MeV LDM decays and for f^bs = 1- The lines are the same as in the left panel. 



by the differential brightness temperature between a neutral 
hydrogen patch and the CMB: 



am Ts — Tcmb 

oTb ~ — r, 

1+z 



(18) 



where r is the optical depth of the neutral IGM at 21(1 -|- 
z) cm: 



327:kBi^o''TsH{z) 



Mi 



(19) 



In equation (19), hp and fcs are the Planck and Boltzmann 
constants, respectively, uq = 1420 MHz is the 21 cm hyper- 
fine transition frequency, and A/hi is the local HI number 
density. If Ts is higher than Tcmb, the neutral IGM will 
be visible in emission against the CMB; on the contrary, if 
Ts < Tcmb it will be visible in absorption. 



4 RESULTS 

In a Universe where DM does not decay or annihilate the 

spin temperature and the kinetic temperature of the gas 
track the CMB temperature down to z ~ 300, when the 
kinetic temperature starts to decrease adiabatically, Tk oc 
(1 + z^, while Tcmb oc {1 + z). The collisions at this red- 
shift are efficient at coupling Ts and Tk'- the spin temper- 
ature subsequently tracks the kinetic temperature down to 
2: ~ 70. At lower rodshifts radiative coupling to the CMB 
becomes dominant and Ts Tcmb again. As a result, for 
30 < z < 300 HI is visible in absorption against the CMB at 
wavelength of 21 (I-I-2:) cm. This scenario could change con- 
siderably if wc allow for DM decay/annihilation, depressing 
or even erasing the absorption feature discussed above. Also, 
we note that different DM candidates leave different traces 
on the 21 cm background signal, and therefore it could be 



possible, in principle, to constrain directly the DM nature 
through 21 cm observations. To isolate the effect of DM, 
we will assume in the following that the entire cosmic dark 

matter content is constituted by particles of the considered 
type (sterile neutrinos or LDM), i.e. ^Ix = — ^b- 



4.1 WDM: 25-keV sterile neutrino decay 

In Fig. 1 we compare the effects of 25-kcV sterile neutrino 
decays on Ts assuming either the physically motivated value 
of fabs (left panel) or the commonly adopted value fabs = 0.5 
(right panel) , i.e. the maximum allowed value for sterile neu- 
trinos corresponding to complete absorption. In addition, in 
each of the panels we explore the differences between mod- 
els with or without the energy injection term due to DM 
decays. We will refer in the following to the latter as the 
standard case, for brevity. 

With respect to the case without DM energy iujection, 
we find, independently of the assumption made for fabs, a 
higher kinetic temperature and hydrogen ionization fraction. 
However, such an enhancement produces only very modest 
Ts differences with respect to the standard case, which has 
some effect on the differential brightness temperature STb. 
The difference Ts — Tcmb remains extremely small due to 
the limited ability of the Lya pumping to decouple the two 
temperatures. 

The main difference between the results obtained from 
the two fabs assumptions is the behavior of Tk at z ~ 10—70. 
When we use the realistic /ats value, the IGM kinetic temper- 
ature positively deviates at most by 50% from its standard 
evolution; if instead fats = 0.5 is adopted (e.g. as assumed 
for example by Furlanetto et al. 2006), the kinetic tempera- 
ture overshoots the CMB one and reaches very large values 
(> 100 K) at low redshifts. 
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Figure 3. Left panel: Ts and Tj, as a function of redshift (10-MeV LDM annihilations, as in RMF06a). Thin solid line: Tcmb', thick 
solid (short-dashed) line: Tg with (without) 10-MeV LDM annihilations; dotted (long-dashed) line: T^. with (without) LDM annihilations. 
Right: Ts and as a function of redshift for 10-MeV LDM annihilations and for f^hg = 1. The lines are the same as in the left panel. 



4.2 LDM: 10-MeV decay/annihilation 

We now turn to the analysis of LDM candidates. Following 
the same procedure as in the previous Section we compare 
the results of the two different assumptions for fabs (note 
that for LDM, the constant absorbed fraction case corre- 
sponds to fabs = 1; see Sec. 2.1), and also test the effects of 
decays/annihilation against the standard case. 

We start from the analysis of the realistic fabs case (Fig. 
2, left panel) for decaying LDM. The most striking feature 
is that the evolution of Tk decouples from the adiabatic one 
already at a « 50, and starts to increase below z « 20, 
reaching Tk = 30 K at a = 10. Such thermal history forces 
Ts to remain below Tcmb for a much longer redshift interval 
(down to z = 10) than in the standard case, in which Ts ~ 
Tcmb already at 2; « 25. Interestingly, this effect extends 
the frequency range in which the IGM can be observed in 
absorption to higher values. The main physical reason for the 
larger impact of LDM decays on Ts with respect to sterile 
neutrinos analyzed above basically lies in the larger fabs 
value (approximately a factor of 10 below z = 30) of LDM 
particles; this is on top of their larger rest mass. The higher 
heating rate increases both Tk and the Lya background {i. e. 
ya), thus pushing Ts away from Tcmb- 

When fabs is instead artificially forced to be constant 
and equal to unity (Fig. 2, right panel), the physical argu- 
ments given above can still be used to interpret the results, 
but the deviation from the standard case is much more dra- 
matic, and results in an overestimate of the LDM decays im- 
pact. For example, Tk becomes larger than Tcmb already 
at 2 « 50 and stays above it thereafter, reaching > 1000 K 
at 2 = 10. The increase is so strong that it essentially erases 
the absorption feature expected in the standard case above 
z = 30. The IGM according to this prescription should be 
observed in emission up to very high redshifts. 



For comparison sake, we note that this case has nearly 
the same decay rate, Tx = 2.5 x 10~^® s~^, as that assumed 
in one of the cases explored by Furlanetto et al. (2006) (dot- 
dashed line in Fig. 1 of that paper). Consistently, the evo- 
lution of Tk from the two studies is essentially the same. 

Finally we turn to the case of annihilating LDM (Fig. 
3). As usual, the left panel of that figure reports the case 
in which the realistic, redshift dependent, values of fabs are 
assumed. In this case, the kinetic temperature deviates from 
the standard adiabatic evolution at extremly high redshifts 
{z > 200): this is because the annihilation process depends 
on the square of the baryon density. As a result, its contribu- 
tion predominantly occurs at early epochs and progressively 
vanishes with time, as realized also from the fact that be- 
low z = 100 the Tk curve is simply shifted to an higher 
adiabat. This fact causes the spin temperature to remain 
closer to (although always lower than) Tcmb, thus preserv- 
ing the standard absorption feature. Quite remarkably the 
results assuming the realistic fabs values do not differ ap- 
preciably from those obtained by imposing fabs = 1 (Fig. 
3, right panel). As pointed out before, the effects of anni- 
hilating DM are evident only at very high redshift, where 
the detailed calculation gives fabs ~ 1 (see RFM06a), thus 
making the non-evolving fabs approximation acceptable. 



4.3 Global 21 cm background. 

Having obtained the evolution of Ts for the different DM 
candidates, we are now ready to compute for each case the 
quantity that is most readily associated with observations. 
I.e. the differential brightness temperature, 5Tb. These are 
shown in Fig. 4 for the realistic (left panel) and the con- 
stant fabs (right panel) cases, respectively. In addition, to 
facilitate the comparison, we have also plotted in Fig. 5 the 
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Figure 4. Left panel: 21 cm differential brightness temperature as a function of redshift. The solid line shows (5T(, without decay- 
ing/annihilating DM; while the long dashed, short dashed and dotted lines refer to STi, with 25-keV decaying WDM, 10-MeV decaying 
LDM and 10-MeV annihilating LDM, respectively. Calculations were performed assuming f^bs as in RMF06a. Right: 21 cm differential 
brightness temperature as a function of redshift as in the left panel. Calculations were performed assuming f^^s = 0.5 for sterile neutrinos 

and faba = 1 in the other cases. 



brightness temperature deviation STb — STb^ of the various 
models from the standard one, restricted to the reaUstic fabs 
evolution only. 

As usual, we start our analysis from the realistic fabs 
case. For sterile neutrinos (long-dashed curves) the charac- 
teristic absorption feature in the 21 cm background signal 
expected at 30 < z < 300 is only slightly modified from the 
standard case (soUd line): the maximum difference is found 
to be only w —2 mK in the range z ~10— 40. Yet, such a 
small signature of WDM decays could still be in principle ob- 
servable: even modest sized single-dish radio telescopes can 
reach the required mK sensitivity in an all sky observation. 
The real challenge for the observation is represented by the 
ability to disentangle this cosmological signal from the vari- 
ous foregrounds (particularly the Galactic synchrotron emis- 
sion) which could be several orders of magnitude brighter. 
The case of decaying 10-MeV LDM (short-dashed fines in 
Fig.4) is more interesting, as the difference with the stan- 
dard case are larger. For this case we find that the deviation 
is predicted to be ~ —(5—8) mK in the range 20 < z < 40. 
Such an amplitude could be detected by LOFAR or SKA 
after a 1000 hour all sky integration and by most single 
dish radio telescopes, provide foreground contamination is 
taken care of. Finally, when LDM annihilations are consid- 
ered, a considerably different 5Tb evolution is obtained with 
respect to the standard case. Larger deviations are present 
in the entire range 40 < z < 200, where the energy in- 
put of LDM annihilations forces 5Tb to values larger than 
—20 mK, about two times smaller (in absolute value) than 
for the standard evolution. Such a large difference could in 
principle facilitate discriminating the annihilation scenario 
from the standard one. In practice, though, foregrounds and 



ionospheric contamination become more severe at the low 
observing frequencies implied by these high redshifts. 

Although of academic interest only, it is instructive to 
compare the constant fabs cases (Fig. 4, right panel) with 
the previous results. As seen clearly form the Figure, con- 
siderably different conclusions would be drawn. First, sterile 
neutrino would drive 5Tb to positive values, resulting in an 
emission signal below z = 20; a similar trend is also obtained 
for LDM decays, whose effects are seen in 21 cm emission 
with an amplitude of about 10 mK. The results for LDM 
annihilation instead do not show appreciable dependencies 
on the fabs prescription, for the reasons discussed above. 



5 DISCUSSION 

From the results obtained in this paper, we conclude that 
it is in principle possible to observe the HI 21 cm signal 
from the Dark Ages produced by the energy input due to 
decays/annihilations of the most popular light/warm DM 
candidates. If so, radio observations might represent one of 
the most promising tools to study the nature of DM, as dif- 
ferent particles are predicted to leave a specific signature 
on the signal. The sensitivity required to measure the 21 cm 
background signal can be achieved not only by the next gen- 
eration of radio interferometers, but also by existing radio 
observatories. 

However, the various foregrounds {i.e. Galactic free- 
free and synchrotron emission, unresolved extra-galactic ra- 
dio sources, free-free emission from ionizing sources, syn- 
chrotron emission from cluster radio halos and relics) are 
much stronger than the cosmological signal, and will cer- 
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tainly prove extremely difficult to remove. Hence, a clear 
detection could be challenging to achieve. 

It is beyond the scope of this paper to deal in detail 
with the problems of ionospheric scintillation and foreground 
contamination. A number of studies have discussed the fore- 
grounds complications in some detail (e.g. Shaver et al. 1999; 
Oh & Mack 2003; Di Matteo, Ciardi & Miniati 2004) and 
we refer the reader to those papers for more information. 
For our aims it is sufficient to say that in general the sky 
temperature can be roughly described as: 



t sky 



180 



USOMHzJ ^ 



(20) 



(see e.g. Furlanetto, Oh & Briggs 2006) and is, at the fre- 
quencies of interest, several orders of magnitude higher than 
the cosmological 21 cm background signal. 

The success of the 21 cm background observations 
to constrain DM will depend on the capability of effec- 
tively removing the foregrounds and of correcting for iono- 
spheric variations (e.g. Gnedin & Shaver 2004; Zaldaxriaga, 
Furlanetto & Hernquist 2004; Santos, Cooray & Knox 2005; 
Morales, Bowman & Hewitt 2005). 

An alternative method to study the 21 cm backgrounds 
is by separating its spectral features from the smooth fore- 
ground spectrum (e.g. Furlanetto 2006, Shaver et al. 1999). 
Fig. 6 shows the gradient of 5Tb for the DM candidates 
studied here. It is clear from the figure that the gradients 
are higher and easier to constrain in the frequency range 
v = 1420/(1 -I- «)MHz ~ 10-40 MHz. The foreground gra- 
dients however increase with decreasing frequency making 
it difficult to predict the best frequencies to get to for a 
successful detection. 

Discriminating among different DM cases would imply 
the ability to distinguish the difference in brightness temper- 
ature and/or in its gradient with respect to the standard sce- 
nario. We define these two quantities as ASTt = STb — 5Tb^ 
and ASTb/Af = dSTb/du-dSTb" /diy, where the stands for 
the standard {i.e. non decaying/annihilating) DM scenario. 

Fig. 7 shows ASTb versus ASTb/Av for the three DM 
models: the blue triangles correspond to 10-MeV annihila- 
tions, the red squares to 10-McV decays and the black stars 
to the 25-keV sterile neutrino decays. The points are sepa- 
rated by 1 MHz steps in a range that goes from 5 to 130 
MHz. As the sensitivity is somewhat uncertain for the dif- 
ferent proposed 21 cm experiments, we assume as a guide- 
line that a successful detection requires ASTb > 3 mK and 
ASTb/Av > 0.6 mK/MHz, respectively. The adopted value 
for ASTb is a conservative one given that mK sensitivities 
are already achievable in all sky observations even by ex- 
isting single dish radio telescopes. It has been chosen so to 
partially take into consideration the difficulties associated 
with foregrounds removal. 

The shaded area covers the parameter space corre- 
sponding to likely non-detections. For the 25-keV sterile neu- 
trinos case the points lie almost entirely within the shaded 
area, thus future observations would hardly be able to pin 
down this scenario. 

In the case of decaying or annihilating 10-MeV LDM 
positive detections could be achieved on different frequency 
ranges. For the decays the ideal frequencies to study would 
be 40-80 MHz, where ASTb ^ 4 mK, and 80-90 MHz, 
where ASTb ~ 2 mK but A<5Tb/Az/ > 0.6 mK/MHz. At 




Figure 5. 21 cm differential brightness temperature difference 
in mK between the three DM cases considered and a non 
decaying— annihilating DM scenario. The dashed, solid and dot- 
ted lines correspond to 25-kev sterile neutrino decays, 10-Mev 
LDM decays and 10-MeV annihilations case respectively. 



1/ > 90 MHz, or 2 < 15 the first luminous sources start 
having a dominant impact on the IGM and the small mod- 
ifications from decaying/annihilating DM would not be vis- 
ible. 

For the 10-MeV LDM annihilations case the best ob- 
serving froquoncios seem to bo 10—30 MHz. At these fre- 
quencies ASTb ^ 6 mK with a peak of over 20 mK at 
10—20 MHz. Observations at such low frequencies arc par- 
ticularly difficult and probably future radio interferometers 
will not reach such large wavelengths. Probably 30 MHz is 
the minimum frequency which will be achieved by next gen- 
eration instruments, and the frequency at which it will then 
be possible to constrain 10-MeV annihilating LDM. 

We conclude then that it is likely that future observa- 
tions of the 21 cm background radiation will allow us to ef- 
fectively constrain among some of the DM candidates, which 
by decays or annihilations can influence deeply the proper- 
ties of the IGM. 

Nevertheless our study shows that accounting for the 
physical evolution of fabs (z) strongly reduces the amplitude 
of DM effects on the IGM with respect to those predicted by 
previous studies: for example we find that the effect of 25- 
keV decaying sterile neutrinos would be hardly observable. 



REFERENCES 

Abazajian K., 2006, PhRvD, 73f, 3506 

Abazajian K., Fuller G. M., Tucker W. H., 2001, ApJ, 562, 593 
Abazajian K., Koushiappas S. M., 2006, Phys.Rev., D74, 023527 
Allison, A. C, Dalgarno, A., 1969, ApJ, 158, 423 
Beacom J. F., Bell N. F., Bertone G., 2005, PhRvL, 94q, 1301 
Beacom J. F., Yuksel H., 2006, PhRvL, 97g, 1102 



© 0000 RAS, MNRAS 000, 000-000 



Constraining DM through 21 cm observations 9 




^6l:J^v (mK/MHz) 



Figure 6. Gradient of tlie 21 cm differential brightness tempera- 
ture. The soUd hne shows dST^^/dv without decaying/annihilating 
DM; while the long dashed, short dashed and dotted lines refer 
to d&T^y/dv with 25-keV decaying WDM, 10-MeV decaying LDM 
and 10-MeV annihilating LDM, respectively. Calculations were 
performed assuming f^hs ^ in RMF06a. 

Biermann P. L., Kusenko A., 2006, PhRvL, 96L, 1301 
Bowman, J., Morales, M., Hewitt, J., 2005, astro-ph/0507357 
Boyarsky A., Neronov A., Ruchayskiy O., Shaposhnikov M., 

2006a, MNRAS, 370, 213, astro-ph/0512509 
Boyarsky A., Neronov A., Ruchayskiy O., Shaposhnikov M., 

2006b, astro-ph/0603368 
Chen, X. & Kamionkowski, M., 2004, PhRvD, 70, 043502 
Chen, X. & Miralda-Escude, J. 2004, ApJ, 602, 1 
Ciardi, B. & Madau, P. 2003, ApJ, 596, 1 

Di Matteo, T., Ciardi, B., Miniati, F., 2004, MNRAS, 355, 1053 

Dolgov A. D., 2002, Phys.Rept., 370, 333 

Field, G. B. 1959, ApJ, 129, 551 

Furlanetto S. R., 2006, MNRAS, 371, 867 

Furlanetto S. R., Oh S. P., Briggs, 2006, PhR, 433, 181F 

Furlanetto S. R., Oh S. P., Pierpaoli E., 2006, PhRvD, 74j, 3502 

Gnedin, N. Y., Shaver, P. A., 2004, ApJ, 608, 611 

Hansen S. H., Haiman Z., 2004, ApJ, 600, 26 

Hirata, CM. 2005, ApJ, in press (astro-ph/0507102) 

Kassim, N. E., Lazio, T. J. W., Ray, P. S., Crane, P. C, Hicks, B. 

C, Stewart, K. P., Cohen, A. S., & Lane, W. M. 2004, Planet. 

Space Sci., 52, 1343 
Knodlseder J. et al., 2005, A&A, 441, 513 
Liszt, H., 2001, A& A, 371, 698L 

Madau, P., Meiksin, A., & Rees, M. J. 1997, ApJ, 475, 492 
Mapelli M., Ferrara A., 2005, MNRAS, 364, 2 
Mapelli M., Ferrara A. & Pierpaoli E., 2006, MNRAS, 369, 1719 
Morales, M. F., Bowman, J. D., Hewitt, J. N., 2005, astro- 
ph/0510027 

Oh, S. P., Mack, K. J., 2003, MNRAS, 346, 871 
Padmanabhan N., Finkbeiner D. P., 2005, PhRvD, 72b, 3508 
Pen, U. L.,Wu, X. P., Peterson, J. 2004, astro-ph/0404083 
Pengelly, R. M., 1964, MNRAS, 127, 145 
Pierpaoli E., 2004, PhRvL, 92, 031301 

Ripamonti E., Mapelli M., Ferrara A., 2006a, MNRAS, accepted. 



Figure 7. Deviations of the differential brightness temperature 
(j/-axis) and its gradient (x-axis) from those predicted in the stan- 
dard case. The diagram illustrates the capability of constraining 
DM with future 21 cm observations. The black stars, red squares 
and blue triangles represent 25-keV decaying WDM, 10-MeV de- 
caying LDM and 10-MeV annihilating LDM, respectively. The 
shaded area indicates the parameter space in which detection is 
excluded given the expected observational capabilities. Numbers 
along the curves refer to frequencies in MHz units. 

astro-ph/0606482 [RMF06a] 
Ripamonti E., MapeUi M., Ferrara A., 2006b, MNRAS, submit- 
ted, astro-ph/0606483 [RMF06b] 
Santos, M. G., Cooray, A., Knox, L., 2005, ApJ, 625, 575 
Seager, S., Sasselov, D. D., Scott, D., 1999, ApJ, 523, 1 
Shaver, P., Windhorst, R., Madau, P., de Bruyn, G., 1999, A&A, 
345, 380 

Shchekinov Y.A., Vasiliev E.O., 2006, MNRAS, submitted, astro- 

ph/0604231 
ShuU J. M., 1979, ApJ, 234, 761 

ShuU J. M., van Steenberg M. E., 1985, ApJ, 298, 268 
Spergel, D. N. et al., 2006, ApJ, submitted 

Watson C. R., Beacom J. F., Yiiksel H., Walker T. P., 2006, 

PhRvD, 74c, 3009 
Wouthuysen, S. A. 1952, AJ, 57, 31 

Wyithe, J.S., Loeb, A., Barnes, D., 2005, astro-ph/0506045 
Zaldarriaga, M., Furlanetto, S. R., Hernquist, L., 2004, ApJ, 608, 
622 

Zhang L., Chen X., Lei Y.-A., Si Z., 2006, PhRvD, submitted, 
astro-ph/0603425 



© 0000 RAS, MNRAS 000, 000-000 



